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Various recent experiments hint at a geometry dependence 
of scaling relations in Rayleigh-Benard convection. Aspect 
ratio and shape dependences have been found. In this paper 
a mechanism is offered which can account for such depen- 
dences. It is based on Prandtl's theory for laminar boundary 
layers and on the conservation of volume flux of the large 
scale wind. The mechanism implies the possibility of differ- 
ent thicknesses of the kinetic boundary layers at the sidewalls 
and the top/bottom plates, just as experimentally found by 
Qiu and Xia (Phys. Rev. E58, 486 (1998)), and also different 
i?a-scaling of the wind measured over the plates and at the 
sidewalls. In the second part of the paper a scaling argument 
for the velocity and temperature fluctuations in the bulk is 
developed. 



Turbulent Rayleigh-Benard convection is one of the 
classical problems in fluid dynamics |l] ^ . The great in- 
terest into this problem presumably also originates from 
the relevance of thermal turbulence in meteorology, geo- 
physics, oceanography, and astrophysics. However, in 
contrast to these natural problems, in laboratory experi- 
ments the thermal convection is confined to a container. 
Recent experiments have revealed that geometrical de- 
tails of this container are important. Daya and Ecke Q 
find that the bulk temperature and velocity fluctuations 
in a cylindrical cell and in a square cell are very different 
from each other. For cylindrical cells the aspect ratio de- 
pendence has been examined by the Ahlers group |^|,|^] 
and by the Tong group 0; in those experiments only 
the prefactors of effective power laws seem to change. 
Niemela and Sreenivasan H] find an aspect ratio depen- 
dence of the heat flow at large Rayleigh numbers, and 
conclude that the focus of the next generation of experi- 
ments must be on large aspect ratio cells. 

Another recent surprising result to be mentioned in 
this context is the dependence of the kinetic boundary 
layer (BL) thickness (which is set by the mean flow ve- 
locity profile above the wall) on its location in the cell. 
Qiu and Xia p|, [I0| find that it scales differently at the 
sidewalls and at the top or bottom plates. They find 
8 W /L = 3.6 i? a -o-26±o.o3 for the si dewall BL thickness 
(from now on abbreviated with w, for wall) and S p /L = 
0.65-Ra-°- 16±0 - 02 Pr - 24±0 - 01 for the thickness at the top 
or bottom plates (from now on abbreviated with p, for 
plate), both for an aspect ratio r — 1 cylindrical cell. No 
existing theory of Rayleigh-Benard convection accounts 
for this difference. Our own unifying theory of ther- 
mal convection [[TTl[L2j so far allows only for one kinetic 



boundary layer thickness, and gives a power law expo- 
nent —0.22 (for r = 1 and a Prandtl number of Pr = 6) 
H, cl ose to the average (-0.16 - 0.26)/2 = -0.21 of 
the experimental results for the plate and the wall BLs 
thicknesses. 

In this paper we set out to offer a mechanism to 
account for the shape and aspect ratio dependence in 
Rayleigh-Benard convection, and for the difference be- 
tween the wall and plate kinetic BLs. The starting point 
is Prandtl's BL theory []R3|-|l5|| . The main physical ingre- 
dient is the conservation of the volume flux of the large 
scale wind. In the second part of the paper we present 
scaling arguments for the fluctuations of the velocity and 
the temperature in the bulk. 

Laminar BL flow is described by the famous Prandtl 
equations ]l3|-|l5[] . In these equations streamwise lengths 
are scaled by the streamwise length scale I, while wall- 
normal lengths are (re)scaled by l/y/Re. Here, Re — 
IUq/u is the Reynolds number based on the streamwise 
length scale / and the streamwise velocity scale Uq. Cor- 
respondingly, streamwise velocities are scaled by Uq and 
wall-normal velocities by Uo/VRe. The immediate con- 
sequence is that the thickness of the laminar BL scales 
as S ~ l/y/Re ( (r|; see also e.g. section 39 of ref. ft5f). 

In the context of Rayleigh-Benard convection this 
means that the relevant streamwise length scales are the 
width d of the cell for the plates and the height L of the 
cell for the walls, i.e., different for cells with non-unity 
aspect ratio r = d/L. Correspondingly, according to the 
Prandtl theory, also the widths of the BLs at the walls 
and at the plates are different, namely, 
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respectively. Here, a is a dimensionless prefactor of order 
1; Blasius fl4[| gave a — 1.72 for a semi-infinite kinetic 
BL, cf. also ]15|. If we assumed that the streamwise wall 
and plate velocities are the same, U w = U v , we would 
immediately get aspect ratio dependent BL widths S w — 

Sp/Vr. 

However, the experiments by Tong's group show 
that the assumption U w = U p does not hold. Indeed, a 
more realistic assumption seems to us volume flux con- 
servation, i.e., equal volume fluxes of the wind over the 
walls and the plates, 
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V W = V P . (3) 

Then the geometry of the cell immediately enters. We 
will discuss two extreme situations: The "confined" flow, 
and the "unconfined" flow, see figure |l|. In the confined 
flow the convection roll has the same spanwise exten- 
sion both at the plates and the walls. In contrast, the 
unconfined flow makes use of the full plate width d in 
the middle, whereas towards the walls its extension is of 
order 2\J d5 w fl6|j . The real flow will presumably be in 
between these two extremes. 

Confined flow: The spanwise length scale of the flow 
is of order 2\/d6 w at the walls and at the plates (see 
figure |l|). The corresponding volume fluxes are V w = 
2^/d5 w S w U w at the walls and V p = 2y/dS w 5J[J p at the 
plates. With volume flux conservation eq. (f|) and the 
Prandtl relations (|l|) and (|^) we immediately get 

s w /s p = r- 1 , (4) 
u p /u w = r- 1 . (5) 

Thus the wall and plate BLs have different thicknesses. 
But there is no difference in the Ra-scaling of S w and 5 p , 
in contrast to what was found by Xia's group j9|. 

Unconfined flow: The spanwise length scale of the flow 
at the cylindrical walls still is 2\/dS w . But now by defi- 
nition the flow at the plates is unconfined spanwise and 
the respective spanwise extension is of order d. The cor- 
responding volume fluxes are now V w = 2y / dS w 5 w Urn at 
the wall and V p = d5 p U p at the plate. With volume flux 
conservation eq. (p|) and the Prandtl relations (|l|) and (|^) 
we now obtain U p /U w = isL~ 1 r~ 4 2 i a A and 

S W /5 P = r-Wy/ujju*, . (6) 

We now introduce wall and plate Reynolds numbers, 

Re w = U w L/v , Re p = XJ p Ljv . (7) 

Note that both are defined with the height L of the cell. 
With these definitions we find 

6 W /6 P = T-^HaRep- 1 / 2 = V^T^Re- 1 ^, (8) 

U p /U w = Re p /Re w = T^l^Re' 1 = AaT' 2 Re~ 1/2 . 

(9) 

Equation (||) is our first main result. It shows, for 
unconfined flow, indeed a Reynolds number and thus a 
Rayleigh number dependence of the ratio between the BL 
thicknesses at the wall and at the plate, due to the cylin- 
drical shape of the container. S w and 5 P scale differently 
with Ra, just as experimentally observed Q. 

Since experimental flow is expected to lie in between 
the confined and the unconfined cases, 5 W /5 P should be 




FIG. 1. Confined (a) and unconfined (b) flow in a cylindri- 
cal cell. The upper drawings show sketches of the 3D flow, 
the lower ones a top view on the bottom plate. The flow area 
above the plate is shaded. The spanwise extension of the wind 
along the walls is 2V dS w and follows from trigonometry. 

between Re p X ^ 2 and i?e° or between Ra~ 0/22 and Ra°; 
here we have used the effective scaling Re p ~ Ra QA5 
which holds for Pr = 5 and Y = I both experimentally 
and theoretically |l2| in the Ra regime I0 8 < Ra < 
10 . Indeed, the experimental ratio has a scaling expo- 
nent in between the values given by the unconfined and 
confined cases, namely S w /5 p ~ Ra~ 011 fe|. 

We now come to our other result, eq.(0), and to the 
aspect ratio dependence of S w /S p and U p /U w . For T = 
1/2 and T = 2 these ratios (as predicted in the two cases, 
and gauged to T = 1) are summarized in table §. The 
experimental data are obtained from figure 3 (r = 1, 
Ra = 3.7 • 10 9 ), figure 6 (r = 2, Ra = 4.9 ■ 10 s ), and 
figure 8 (T = 1/2, Ra = 3.28 • 10 10 ) of reference @. 
In this series of figures both T and Ra were changed at 
the same time. In order to compare the ratios at T = 2 
measured at Ra — 5- 10 8 with those at T = 1 measured at 
Ra = 3.7 • 10 9 , they first have to be extrapolated to this 
higher Ra number. Correspondingly, in the T = 1/2 case 
we have to extrapolate the ratios from the original value 
Ra = 3.3 • 10 10 to Ra = 3.7 • 10 9 . These extrapolations 
have been done according to the unconfined case formulas 
(||), (||) or according to the confined case ones (|J), (ph. 
In the latter case in fact no extrapolation is necessary as 
in those equations there is no Ra or Re dependence. 

From table || we conclude that the experimental re- 
sults (if extrapolated with the unconfined flow model) 
in general lie in between the predictions of the confined 
and unconfined cases, just as expected. An exception is 
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model 


r = 1 


r = 1/2 


r = 2 


Sw/Sp 


experi., extrapol: confined 


1 


2.1 


0.29 




experi., extrapol: unconfined 


1 


3.3 


0.19 




theory, confined 


1 


2.0 


0.50 




theory, unconfined 


1 


5.6 


0.18 


u p /u w 


experi., extrapol: confined 


1 


1.2 


2.0 




experi., extrapol: unconfined 


1 


3.1 


0.8 




theory, confined 


1 


2.0 


0.5 




theory, unconfined 


1 


16 


0.063 



TABLE I. S w /S p and U p /U w = Re p /Re w for different as- 
pect ratios, gauged to T = 1, where we have put these ratios 
to 1. The Reynolds number Re p = 3.7 ■ 10 9 is fixed and the 
experimental values are extrapolated to that Re p , see text. 
Note that the error in reading off in particular the bound- 
ary layer widths from the figures 3,6, and 8 of ref. j?| is at 
least 25% and therefore the experimental ratios should only 
be taken as a very rough estimate. 



U p /U w for r = 2; we have no explanation for that. If 
extrapolated with the confined flow model both velocity 
ratios do not lie in between the predictions of the two 
flow cases. 

The physics of the aspect ratio dependence is as fol- 
lows: For increasing aspect ratio the available plate area 
of the flow increases and the plate velocity can go down. 
According to the Prandtl law (||) , the plate BL thickness 
then increases. At the walls there is no additional lateral 
space for larger aspect ratios, therefore the flow has to 
accelerate, U w > U p , and correspondingly S w < S p . 

Note that above dependence on T was derived assum- 
ing there is one convection roll. Once the convection 
develops two rolls on top of each other (for T less than 
1) or next to each other (for V larger than 1), the aspect 
ratio dependence has to be embodied in a more sophisti- 
cated way. For small T < 1/2 a transition from one roll 
to two rolls has been observed numerically pi| , [l9[] . Pos- 
sibly also in experiment both states of convection have 
been realized, see the discussion in || . Thus one cannot 
formally consider T — > or T — ► oo in the relations (Q) - 

Can one manipulate a system such that it is no longer 
in between the confined and the unconfined cases? An 
option may be to take a square box and slightly tilt it 
in order to disfavor the flow along the diagonal, which 
would occur without tilting Q. If this can be achieved, 
one expects such flow to follow the confined flow model 
with confinement by the left and right walls. For this type 
of flow eq. (Q) should hold, i.e., the wall and the plate BLs 
should show the same i?a-scaling. An experimental test 
of this predictions seems worthwhile to us to check the 
proposed theory. 

Kinetic and thermal dissipation in the BLs: Within 
above framework the aspect ratio dependence can also 
be embodied into the unifying theory of thermal convec- 



tion [pdl,|l2|| . The main idea of that theory is to split 
the kinetic dissipation e u and the thermal dissipation eg, 
for which exact global relations can be derived from the 
Boussinesq equations, into their bulk and BL contribu- 
tions. 

For the kinetic dissipation this splitting obviously has 
to be extended to 



+ e«,6 



(10) 



where e UlP , e UiW mean the kinetic dissipation in the plates' 
or walls' BLs, and e u f, the bulk kinetic dissipation. Fol- 
lowing pd[, the BL dissipations are estimated as 



and 
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V 1 P r 5/2 _ ^9 5 o 3 / 2 r~ 9 / 2 /?p 5 / 4 
L 4 aT 1 / 2 p ~ L 4 w 

(11) 
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The bulk dissipation will be dealt with below. 

For the thermal dissipation the partition into BLs and 
bulk is different. Though at least in ideal RB convection 
the sidewalk are perfectly isolated (no flux condition), 
thermal BLs can develop at the sidewalls, as observed 
both in numerical simulations and in experiment [fflJly- 
plj . However, within our unifying theory ref. Tl|.|12[ ^ 
is not necessary to distinguish between sidewall thermal 
dissipation and bulk thermal dissipation. Therefore, here 
we introduce the notation eg tP for the thermal dissipation 
in the plate BLs and eg t p for the thermal dissipation else- 
where, i.e., within the bulk and within the sidewall BLs 
( "non-plate" ) . 

To estimate these contributions we remind that eg = 
kA 2 L~ 2 Nu and that the heat current Nu consists of two 
terms 

Nu = (K0) A)t (z) - nd z {6) A>t (z)) /(kAL- 1 ) 

= {ee,p + ee, P ) /(kA 2 L- 2 ) = e 8 /( K A 2 L- 2 ). (13) 



Here 



I A.t 



denotes the average on time and x, y-plane. 



For z outside of the plates' thermal BLs the first term 
in eq. ( |l3| ) dominates. We estimate (u z 9) A t ~ UA and 
obtain eg } p/ (kA 2 / L 2 ) ~ PrRe. Note again that the ther- 
mal dissipation in the sidewall BLs is included in eg^ p and 
may even be dominant as compared to the contributions 
from the bulk, as the numerical simulations by Verzicco 
and Camussi suggest Jl9| . 

For z within the plates' thermal BLs, i.e., z/L « 
or z/L w 1, the second term in eq. ( jl3| ) domi- 
nates. Estimating the temperature gradient as A/Xg 
gives eg iP / '(kA 2 / ' L 2 ) ~ L/Xg. To connect the BL width 
Xg with Re, we use the temperature equation dtO = 

1 ~ nXg 2 , which im- 
cf. fTl]] . According to 



UjdjO + ndjdj9, leading to Ud 
plies eg, p /(KA 2 /L 2 ) - y/RePr, 
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this derivation the Reynolds number one has to use in 
eg iP is Re w , while eg tP is determined by Re p . 

Scaling relations for bulk velocity and temperature fluc- 
tuations: Up to now we dealt with the BLs. As for 
the bulk fluctuations u' , the experiments (To) suggest 
that they scale differently from the large scale veloc- 
ity U , namely with a weaker Ra dependence. Typi- 
cally, the fluctuation Reynolds number is Re 1 = vl Ljv ~ 
^ a o.4o±o.03 anc j foe large scale velocity Reynolds num- 
ber with Re = UL/v ~ i? a O-43...o.50 gg^HH. In 
refs. [plO] we developed a theory for the dependences 
Re(Ra, Pr) and Nu(Ra, Pr). Here we add how the fluc- 
tuations Re 1 and 9' should behave as functions of Ra and 
Pr. 

The turbulence in the bulk is driven by the large scale 
wind, therefore e u ^ ~ U 3 / L as used in refs. p| , p^| . The 
energy cascades down and will be dissipated on scales 
comparable to the Kolmogorov scale r/ = [y 3 j eu.b) 1 ! 4 . 
Thus we can estimate the bulk dissipation also with 
vu' 2 /rj 2 . Combining both we obtain 



Re' ~ Re 3 ' 4 



(14) 



With Re ~ i? a 0-43...0.50 this meang Re i _ J?a 0.32...0.38 j 

which is close to above quoted experimental finding 
Re' ~ i? a °- 40±0 03 by Lam et al. p| and also to Daya 
and Ecke's finding g Re' ~ ^0-36^05 for a square ce u. 

The same reasoning can be followed to calculate the 
scaling of the typical bulk temperature fluctuations 9'. 
The thermal bulk dissipation is driven by the large scale 
temperature A; its time scale again is L/U. Thus 
ee b ~ A 2 U/L. On the other hand, the thermal en- 
ergy in the bulk is dissipated by the temperature fluc- 
tuations 9' on scale rj e = (k 3 /£„,&) 1/4 ~ LPr- 3 / 4 Re- 3 / 4 , 
i.e., ee,b ~ k9' 2 /rjg. Balancing these two expressions for 
eg t one obtains 



9'/A ~ Pr-^ 4 Re- 1/4 



(15) 



With the experimental result @ Re ~ Ra 0A3 Pr 
which holds in 10 8 < Ra < 3 • W° and 3 < Pr < 1200, 
this means 9'/A ~ Pr-° 06 Ra~ 011 . For Re ~ Ra 1 / 2 
H,|U one gets 6'/ A ~ Ra-° 13 . Both power laws are 
very close to the experimental findings 9' / A ~ Ra^ 1 / 7 
H or 9' I A ~ Ra-° 1Q±0 02 § for cylindrical cells. 

In conclusion, we have offered a possible mechanism to 
cope with the recently observed effects of the RB cell ge- 
ometry and the different Ra behaviors of the kinetic BLs 
near the plates and the sidewalls. This has to be con- 
firmed in more experimental detail. In particular, the 
comparison of RB convection in a square cell with diag- 
onal (unconfined case) and with edge parallel (confined 
case) wind seems promising. - Another aspect is to dis- 
tinguish the 1-roll from the 2-rolls (or more) states. It 
seems also clear that measuring the large scale wind ve- 
locity should take notice that possibly U p ^ U w . Circu- 
lation time measurements average over both, while local 
methods can distinguish U p and U w . - Finally, we have 



indicated how the global theory of Nu versus Ra, Pr has 
to be modified to take care of geometry effects and how 
the different scalings of bulk fluctuations and large scale 
quantities are related. 
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